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ABSTRACT 

Context. The H30 + molecule probes the oxygen chemistry and the ionization rate of dense circumnuclear gas in galaxies. 

Aims. Recent H30 + observations show variations in the cosmic-ray ionization rate by factors of >10 within our Galaxy. 

Methods. Using the JCMT, we have observed the 364 GHz line of p-H 3 + in the centers of M82 and Arp 220. 

Results. In Arp 220, the line profile suggests that the emission originates in the Western nucleus. In M 82, both the eastern molecular 

peak and the circumnuclear region contribute to the emission. The derived column densities, abundances, and H30 + / H2O ratios 

indicate ionization rates similar to or even exceeding that in the Galactic Center. 

Conclusions. Model calculations of the chemistry of irradiated molecular gas indicate a likely origin of this high ionization rate in 
the extended, evolved starburst of M82. In contrast, irradiation by X-rays from the AGN disk is the most likely model for Arp 220. 
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1. Introduction 

Starbursts and active galaxies often host large masses of molec- 
ular gas and dust in their inner kpc. Gravitational interaction or 
merging cause the funneling (via bars) of the gas towards the 
center. The average gas surface- and number densities are sig- 
nificantly higher (by several orders of magnitude) than in the 
more extended disks. The nuclear molecular gas is often found 
in tori, feeding the activity and thus controlling its evolution. The 
triggering and turn-off mechanisms of the activity are dependent 
on the properties of the gas, which are still poorly known, in 
particular the fraction of the gas at high temperature and den- 
sity. Observations of the molecular gas properties are starting to 
show great promise not only in allowing us to model the evolu- 
tion of the starburst, but also in helping to identify and analyze 
the type of activity lurking in very deeply obscured galactic nu- 
clei. Here, dust extinction is high enough to quench emission 
even down to the mid-inf rare d. For reviews of t his su bject, see 
ISanders & Mirabelldl996l ) and lGao & Solomon! d2004l) . 

The star formation rates in the nuclei of nearby galaxies are 
measured best at (sub-)millimeter wavelengths, where extinction 
is negligible and where a rich spectrum of dust continuum and 
molecular lines allows us to measure fundamental parameters of 
star-forming matt er such as mass, temperature, and density (e.g., 
Aalto et al. 2002). However, existing observations are mostly at 
millimeter wavelengths and thus limited to relatively low tem- 
peratures and densities. The excellent transmission of sites like 
Mauna Kea and Chajnantor enables observation of the HsO + 
molecule, which has a twofold importance. 

First, the H2O molecule acts as a natural filter for select- 
ing hot (^100 K) gas, and therefore traces different regions than 
CO and other molecules commonly studied from low-altitude 



sites. The chemistry of this fi ltering is the evaporation of icy 
grain mantles at T * 100K (I Van der Tak et al.ll2006ah . Other 
molecules evaporating from grain mantles are not as abundant 
as H2O. The H30 + molecule is an excellent proxy for H2O, 
which cannot be observed in thermal lines from the ground at 
redshifts below 0.01. In contrast to HDO whose ratio to H2O 
strongly depends on ambient conditions, the chemical relation 
between H30 + and H2O is well-understood an d the relevant 
reacti ons have been measured in the laboratory (Phillip s~et al.l 
1992). Thus, observations of HaO + are useful probes of the gas 
distribution and kinematics in the warm and dense nuclei of ac- 
tive galaxies. 

Second, the high proton affinity of H2O makes H30 + a key 
ion in the oxygen chemistry of dense molecular gas. Combined 
with information on the H2O abundance, observations of H30 + 
may be used to trace the ionization rate of dense circumnuclear 
clouds by cosmic rays (produced in supernovae, i.e., starbursts) 
and X-rays (from an AGN). The ionization rate of this gas is 
a fundamental parameter for its dynamics, because it regulates 
the efficiency of magnetic support against gravitational collapse. 
Recent APEX observations of H30 + show that the cosmic-ray 
ionization rate in the Sgr B2 cloud near the Galactic center is 
~10x higher than in star-form ing regions at 1-4 kpc from the 
Sun dVan der Tak et alj|2006bl) . 

Motivated by the results of H30 + observations in our Galaxy, 
we searched for HsO + emission in two prototypical active galax- 
ies: M82 and Arp 220. At a distance of 4Mpc, M82 is the pro- 
totype evolved starburst, with a large H2 mass, where many 
molecules were seen for the first time outside our Galaxy, in- 
cluding CO + , a tra cer of strong irradiation of molecular gas 
dFuente et al.l l2006). The molecular emission can be roughly di- 
vided into two prominent lobes of molecular emission, separated 
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by 30", which is interpreted as a 400 pc rotating molecular ring. 
There is also molecular gas in the nucleus, which is best seen in 
high- ./ CO lines or in high-den sity tracer emission, such as HCO 
(e.g jGarcfa-Burillo et al.ll2002l) . 

The prototypical ultraluminous galaxy Arp 220 
(d - 72Mpc), where wate r lines have been det ected 
dGonzalez- Alfonso et al] l2004t ICernicharo et all 120061) . is 
a merger of two gas-rich disk galaxies. Significant (; 1O 9 M ) 
molecular gas has gathered in two nuclei (east and west) that 
have not yet merged into one. The two nuclei are separated 
by 1" and are surrounded by a 5" extended disk of molecular 
gas. The bulk of the molecular emission from Arp 2 20 is 
emerging from this region (e.g., ISakamoto et al.l Il999l) . and 
the JCMT beam covers all of it. Arp 220 has a rich chemistry 
with surprisingly bright HNC 3-2 line emission, outshining 
HCN 3-2 by a factor of 2, which indicates unusual chemical or 
excitational circumstances. High-resolution observations reveal 
that the HNC emission is emerging fr om warm gas (7\ > 40 K) 
near the western nucleus (Aalto 2007). 



2. Observations 

The 3J— 2j line of H30 + was observed during April 26-28, 
2007, with the James Clerk Maxwell Telescope (JCMTf] un- 
der program M07AN03. In the 345 GHz window, the telescope 
has a beam size of 14" and a main beam efficiency of 70%, 
as determined by the telescope staff and verified by scans of 
Mars at the end of each night. One of the central pixels of the 
Heterodyne Array Receiver Program (HARP) was used as front 
end, with receiver temperatures of ~120K at our observing fre- 
quency of ~360GHz. The back end was the Auto-Correlation 
Spectrometer and Imaging System (ACSIS), providing 1.0 GHz 
bandwidth in 2048 channels. The weather was excellent through- 
out, with water vapour columns of ; 1 mm, corresponding to a 
zenith opacity of r! 0.05 at 225 GHz. System temperatures were 
300-400 K for M82 and 200-300 K for Arp 220, which reaches 
higher elevations and whose larger redshift puts the line at a 
frequency with better atmospheric transmission. Double beam 
switching at a rate of 1 Hz was used with an offset of 120" in 
azimuth, which led to flat, stable baselines. Integration times 
(on+off) are 8.2 hrs for Arp 220 and 6.5 hrs for M82, resulting in 
rms noise levels of T* A = 2-3 mK per 8 km s _1 channel. Telescope 
pointing was checked every hour on the CO line emission of 
nearby AGB stars and always found to be within 2". The tele- 
scope z-focus was checked on bright AGB stars at the beginning 
of the night (after sunset) and around midnight, and corrected by 
0.1-0.2 mm. The bright AGB stars were also used to take stan- 
dard spectra, and the results lead us to believe that the calibration 
of the data is accurate to ~10%. 



Table 1. Measured line parameters with lcr errors in brackets in 
units of the last decimal. 



Source 



a(J2000) 
hh mm ss 



<5(J2000) 
/ // 



n 

mK 



V he , 
km s 



AV 
km s 



M82 09:55:52.19 +69:40:48.8 15(2) 269(2) 43(5) 
broad component: 6.5(6) 224(11) 261(23) 

Arp 220 15:34:57.23 +23:30:11.5 4.1(6) 5473(8) 113(20) 




200 

Vheiio (km s" 1 ) 




5400 S500 
v heiio (km s" 1 ) 



Fig. 1. Spectra of the H 3 + 3J-2 2 line toward M82 (top) and 
Arp 220 (bottom), observed with the JCMT. 



dates at this frequency, and the attenuation of the image side- 
band is better than lOdB. After smoothing the data to a velocity 
resolution of 8-16km s _1 , line parameters were extracted by fit- 
ting Gaussian profiles to the spectra. The final spectra are shown 
in Fig. Q]and the line parameters listed in Table [1] Throughout 
this paper, velocities are in the heliocentric frame and redshift is 
computed using the radio convention. 



3. Results 



Initial reduction of the data was done with the Starlink pack- 
age, and the final analysis performed in Class. We inspected the 
individual scans, summing the good ones, and subtracted linear 
baselines. The identification of the spectral featur e as H^O + is 
secure, since the standard molecular line catalogs dPickett et alj 
1998; Mul ler et al.ll200l do not list any other plausible candi- 



1 The JCMT is operated by the Joint Astronomy Center on behalf of 
the Science and Technology Facilities Council of the United Kingdom, 
the Netherlands Organization for Scientific Research, and the National 
Research Council of Canada. 



Emission due to Hj,0 + is clearly detected towards the nuclei of 
both M82 and Arp 220. The signal-to-noise ratio on the line 
profiles is high enough to address the location of the emission 
within the galaxies. In the case of Arp 220, th e line profile is 
single -peaked and similar to that of CH3OH (Martin-Pintado 
2007). The emission peak at Vh e i~5500km s _1 indicates an ori- 
gin in the western nucleus, which is also where th e HNC 3- 
2 em ission peaks in recent interferometric images (Aalt o et all 
2007). In the case of M82, the emission peak at Vh e i~300km s _1 
suggests an origin in the eastern molecular peak at about 12" 
(240 pc) from the nucleus, which is where the CO emission 
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peaks dTilanus et all 19911) and where CO + and CH3OH are also 
located dFuente et all 120061: iMartm et al.ll2006l) . The profile al- 
most extends to Vh e i~200km s _1 , which suggests that signifi- 
cant H30 + is also pres ent in the central region of the galaxy, as 
found before for N2H + dMauersbe rger & Henkelll99lh . In fact, a 
single Gaussian profile does not give a good fit to the M82 spec- 
trum. Figure Q] shows a decomposition of the profile into two 
Gaussians, where the narrow component represents H30 + in the 
eastern molecular peak of M82 and the broad component either 
the more widely-distributed molecular gas or the contribution 
from the nucleus. This decomposition of the line profile may not 
be unique, but serves to illustrate the existence of multiple HiO + 
components in the nucleus of M82. Although the JCMT beam is 
too small to cover the eastern lobe of M82 while pointing at the 
center, it is possible that we are picking up emission from the 
inner side of the lobe. The other HARP pixels do not cover the 
lobe either, but future fully-sampled maps will reveal the spatial 
distribution of H 3 + in M82. 

To estimate the column density of H30 + towards our sources 
from the integrated line intensities (Table [U, we adopt a main 
beam efficiency of 70% and regard the expected kinetic tem- 
perature of ~100K as an upper limit to the excitation tempera- 
ture (r ex ). For r ex = 50 - 100 K, the column density of p-H 3 + 
is 0.7 - 1.3x10° cnr 2 for Arp 220 and 3.0 - 5.5xl0 13 cnr 2 
for M82. The total beam-averaged column density of H30 + is 
2-3 times higher, because the ortho to para ratio of H30 + 
drops from the high-temperature limit of o/p = 1 at 100K 
to o/pl2 at J;50K. The following di scussion adop t s the inter- 
mediate case of o / p — 1.5 ; we r efer to Phil lips et all (1 19921) and 
IVan der Tak et al.l (l2006bl) for further discussion of H30 + radia- 
tive transfer. 

Unlike for lower-frequency lines of heavier molecules 
that are excited by collisions, the excitation temperature of 
H30 + probably reflects the colour temperature (7c) of the 
far-infrared radiation field rather than the kinetic temperature 
of the molecular gas. The low reduced mass of the H30 + 
molecule mak es its excitation very s ensitive to radiative pump - 
ing by dust dPhiUips et alj 119921: IVan der Tak et all |2006b). 
Calcul ations with the RADEX non-LTE radiative transfer pro- 
gram ((Van der T ak et al ] 120071) u sing m olecular data from the 
LAMDA database dSchoier et aT] 120051) indicate that the cou- 
pling between r ex and Tc holds as long as Tq > 30 K and that 
the measured line is optically thin under these conditions. We 
assume a filling factor of the dust radiation field of unity, repre- 
senting the case that the gas and dust are well-mixed. 

The beam-averaged H30 + column density in Arp 220 is al- 
most certainly an underestimate of its source-averaged value, be- 
cause the emission is unlikely to be as extended as 14". First, 
the molecular gas traced by the CO 2-1 emis sion is concen- 
trated in two nuclei of diameter 300 pc each (Sakamot o et al 



Table 2. Source-averaged column densities [cm 2 ] of H2, H2O, 
and H 3 Q + . 



1999). Recent high-resolution work on HNC 3-2 dAalto et al 



2007) further shows that the dense gas is concentrated towards 
the western nucleus with a source size of 0.7 arcsec (260 pc). 
Second, since H30 + is excited by infrared radiation, its emis- 
sion may follow th e infrared light dist ribution. From 3-25 //m 
Keck observations, ISoifer et all (11999) find a size of 0.73" at 
25 /im. They derive a size of 2" at 100 pm, which seems a firm 
upper limit to the H30 + source size. Furthermore, the absence of 
a double-peaked profile from the two nuclei in our data suggests 
that the emission is located on one of the two nuclei and does not 



Source N(H 2 ) 


N(H 2 0) iV(H 3 + ) 


Refs. 


M82 lxlO 23 


<2xl0 14 l.lxlO 14 


[1,2,3] 


Arp 220 lxlO 24 


2xl0 17 l.OxlO 16 


[4,5,6] 



References: [1] Mao et al. (2000); [2] Wilson et al. (2007); 

[3] This work; r41 Downes & Eckart (2007); 
[5] Gonzalez- Alfonso et al. (2004); [6] This work. 



cover the full extent of the CO emission. We thus adopt a source 
size of 0.7", which must of course be verified at high resolution. 

In the case of M82, the dense molecular gas is distributed 
over ma ny clumps as shown by h igh-resolution CO and HCO 
images dGarcia-Burillo et aLll2002l) . These clumps are grouped 
in three lobes (east, middle, and west) separated by «T5" 
(300 pc), with significant emission coming from the medium in 
between. Our JCMT beam is pointed to a position close to the 
central lobe, so that a significant fraction of the beam is probably 
filled with FJ.30 + emission. 

Table [2] reports the column densities of F/_30 + , H2O, and H2 
in our sources, as taken from this work and the literature. The 
values refer to a ^15" area for M82 and a *0.7" area for Arp 
220. Our adopted N(U 2 ) toward Arp 220 is lOx below the peak 
value reported by iDownes & Eckartl d2007l) for a 0.2" area, so 
that it is representative of our assumed H30 + source size of 0.7". 
The table shows that the abundance of H30 + relative to H2 is ~2- 
lOxlO" 9 in these sources. These values are muc h higher than 
the ty pical value for Galactic star-forming cores dPhiUips et al.l 
1 19921) and at or above the value of 3xl0~ 9 in t he Sgr B2 cloud 
at the Galactic center dVan der Tak et a l. 2006b). The ionization 
rate of the molecular gas in the nuclei of Arp 220 and M82 thus 
seems comparable to the value in the Galactic center. This con- 
clusion is supported by the H3OVH2O ratio of 1/50 that we find 
for Arp 220, which is similar to the value of 1/20 for the Galactic 
center. No useful estimate of this ratio can be made for M82: the 
Odin data quoted in Table|2]have a beam size of 2' and thus sam- 
ple a very different gas volume than our JCMT data. Note that 
our estimated H_30 + abundance in Arp 220 is somewhat ham- 
pered by source size uncertainties. 



4. Discussion 

The column densities of H^O + , H2O, and H2 in Table [2] and 
their ratios are compared to chemical models of clouds irra- 
diated by either far-UV or X-ray photons. Far-UV dominated 
clouds are called photon-dominated regions (PDRs) and the 
other type X-ray dominated regions (XDRs). The thermal and 
chemical structures of PDRs and XDRs are quite different, be- 
cause the absorption cross section for X-rays is much lower than 
for far-UV. The ionization fraction in XDRs can be as high as 
x e = 10~ 2 - 10 -1 , while in PDRs, it is 10~ 4 at most. Furthermore, 
in XDRs, larger parts of the cloud can be kept at a high tem- 
perature. For an elaborate dis cussion of the differences b etween 
PDRs and XDRs we refer to J vleiierink & Sp aans (2 0051) where 
the co des are described. The grid of models by (Meiierink et al.l 
2007) is used to look for the best match to our observations. For 
the PDR models, we also consider an elevated cosmic-ray ion- 
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ization rate of £ = 5.0 x 10 15 s '.As refe r ence, we take the 
radiation field found bv lSpaans & Meiierinkl (120071) for M 82 to 
model the CO + abundance found bv lFuente et alj d2006l) . Their 
best-fitting model was an XDR with a density of n = 10 5 cm" 3 
and an incident X-ray flux of Fx — 5.1 erg s" 1 cm" 2 (Go=10 3, ). 
Using this density and radiation field, the model predicts the col- 
umn densities listed in Table [3] The chemical structure of the 
model is shown in Fig. [2] Models for volume densities far be- 
low n = 10 s cm" 3 fail to produce appreciable H30 + abundances, 
which makes it unlikely that the origin of the observed emission 
is in interclum p gas, seen e.g. as the Nh=10 20 cm" 2 CO compo- 
nent in M82 dMao et al.ll2000l) . 

Table [3] shows that the H30 + column density observed in 
M82 may be produced in either XDR or in the high-^ PDR. The 
H3OVH2 ratio in M82 is overproduced in the XDR models and 
underproduced in the standard PDR, but agrees very well with 
the high-£ model. None of the models reproduce the observed 
H3OVH2O ratio in M82, probably because the observations re- 
fer to very different volumes. We conclude that the M82 data are 
matched best by the PDR model with enhanced cosmic-ray flux. 
In contrast, the X-ray models do best in reproducing the obser- 
vations of Arp 220. The H30 + column density of 10 16 cm" 2 is 
matched by an XDR model with /y#=10 23 cm" 2 , scaled up by a 
factor 10, for either value of Fx - The H30 + abundance of 1 x 10" 8 
is matched by the high-fx model for low Nh or the low-Fx 
model for high Nh- The H30 + /H20 ratio is matched by an XDR 
of either Fx, where low values of Nh are preferred. Nevertheless, 
a contribution by the eastern nucleus is implausible because of 
its low mass of dense gas and its low infrared brightness. 

In the context of galactic nuclei, a PDR model represents a 
starburst, where the far-UV is produced by young massive stars 
and the cosmic rays in supernova remnants. The pure PDR thus 
stands for a young starburst and the high-^" model for an evolved 
one. The XDR models represent the much harder radiation envi- 
ronment of an AGN, such as the accretion disk of a supermassive 
black hole. We find that the observations of M82 are matched by 
a high-£ PDR, i.e., an evolved starburst. This conclusion is con- 
sisten t with earlier models of this source dForster Schreiber et alj 
120031) . Our H30 + obs ervations do not require X -ray irradiation 
as in the case of CO + dSpaans & Meiierinkll2007h . 

The fact that the XDR models give conflicting requirements 
on the Nh value toward Arp 220 may well be due to the un- 
usual geometry of this object. The emerging picture of the west- 
ern nucleus of Arp 220 is that there is a 35 pc optically thick 
dust disk surrounding what could potentially be an active nu- 
cleus. The dust disk is warm (To ~ 170 K) and has a peak 
H2 column density of ^10 25 cm" 2 . Hence X-rays, optical, and 
even infrared emission fr om the nucleus itself will be absorbed 
(Downes & Eckart 2007). Surrounding this warm dust disk is a 
cooler (7\ ss 50 K) molecular disk or torus causing CO 2-1 ab- 
sorption features towards the hotter nuclear dust disk in the line 
of sight. Most of the observed starburst activity appears to be as- 
sociated with the rotating disk/torus surrounding the nuclear dust 
disk. Future modeling efforts of H30 + in Arp 220 should take 
this special geometry into account. However, high-resolution ob- 
servations are needed to constrain these models. It is possible for 
example that the disk causes both H30 + emission and absorption 
along the line of sight. 

In the near future, the HIFI heterodyne spectrometer onboard 
ESA's Herschel space observatory will be able to measure H2O 



Table 3. Results of PDR / XDR models. 



N H 

-2 

cm 


-2 -2 

cm cm - 


AT(H 2 0) 
cm" 2 


H 3 + /H 2 


H 3 + /H 2 




XDR F x = 


5.1 erg s" 1 


cm" 2 






lxlO 22 


1.4xl0 21 6.1xl0 13 


5.6xl0 15 


4.3xl0" 8 


1.1x10 


-2 


3x10 


7.9x10 1.9x10 


2.2xl0 16 


2.5xl0" 8 


8.7x10" 


-3 


1 V 1 n 23 

1 A1U 


^ Qv 1 022 1 £ v 1 O 14 
J.7XIU j.DaIU 


5.8xl0 16 


9.2xl0" 9 


6.1x10 


-3 




XDR F x = 


16 erg s" 1 


cm -2 






lxlO 22 


1.7xl0 20 3.1xl0 12 


1.7xl0 14 


1.8xl0" 8 


1.8x10" 


_2 


i w i r\22 

3x10 


2.3x10 1.5x10 


1.2xl0 16 


6.6xl0" 8 


1.2x10" 


-2 


1 v 1 o 23 

1A1U 




8.9xl0 16 


3.4xl0" 8 


9.4x10" 


-3 




PDR G =10 35 


£ = 5.0x lO-'V 






lxlO 22 


4.6X10 2.3X10 


1.6xl0 14 


5.0x10"" 


1.4x10" 


-3 


3xl0 22 


1.5xl0 22 1.7xl0 12 


1.6xl0 15 


1.2xl0" 10 


1.1x10 


-3 


lxlO 23 


5.0xl0 22 6.8xl0 12 


6.6xl0 15 


1.4xl0" 10 


1.0x10" 


-3 




PDR G =10 35 


£ = 5.0x 10- 15 s-' 






lxlO 22 


4.6xl0 21 9.2xl0 12 


2.7xl0 15 


2.0x10"" 


3.4x10" 


-3 


3xl0 22 


1.4xl0 22 4.3xl0 13 


1.3xl0 16 


2.9x10"" 


3.2x10" 


-3 


lxlO 23 


4.9xl0 22 1.6xl0 14 


5.1xl0 16 


3.2x10"' 


3.1x10" 


-3 



and H?0 + lines around 1 THz and make accurate estimates of 
the H20/H?0 + ratios in many sources. These data will allow a 
much better understanding of the ionization rates of molecular 
regions inside and outside our Galaxy. Future additions to the 
model include the effects of freeze-out onto and desorption from 
dust grains, the effect of metallicity changes, and (especially for 
the Arp 220 merger) the time dependence introduced by shocks. 

Acknowledgements. We thank the staff of the JCMT for their support, in 
particular Remo Tilanus, Jan Wouterloot, and Jim Hoge, and thank Rainer 
Mauersberger for the quick and useful referee report. 



References 

Aalto, S. 2007, in Science with ALMA, ed. J. Cernicharo, in press 

Aalto, S., Polatidis, A., Hiittemeister, S., & Curran, S. 2002, A&A, 381, 783 

Aalto, S., Spaans, M., Wiedner, M. C, & Hiittemeister, S. 2007, A&A, 464, 193 

Cernicharo, J., Pardo, J. R., & Weiss, A. 2006, ApJ, 646, L49 

Downes, D. & Eckart, A. 2007, A&A, 468, L57 

Forster Schreiber, N., Genzel, R., Lutz, D., & Sternberg, A. 2003, ApJ, 599, 193 
Fuente, A., Garcfa-Burillo, S., Gerin, M., et al. 2006, ApJ, 641, L105 
Gao, Y. & Solomon, P. M. 2004, ApJ, 606, 271 

Gara'a-Burillo, S., Martin-Pintado, J., Fuente, A., et al. 2002, ApJ, 575, L55 

Gonzalez-Alfonso, E., Smith, H., Fischer, J., et al. 2004, ApJ, 613, 247 

Mao, R. Q., Henkel, C, Schulz, A., et al. 2000, A&A, 358, 433 

Martin, S., Martin-Pintado, J., & Mauersberger, R. 2006, A&A, 450, L13 

Martin-Pintado, J. 2007, in Molecules in Space, ed. F. Combes, in press 

Mauersberger, R. & Henkel, C. 1991, A&A, 245, 457 

Meijerink, R. & Spaans, M. 2005, A&A, 436, 397 

Meijerink, R., Spaans, M., & Israel, F. P. 2007, A&A, 461, 793 

Miiller, H, Thorwirth, S., Roth, D., & Winnewisser, G. 2001, A&A, 370, L49 

Phillips, T. G, van Dishoeck, E. E, & Keene, J. 1992, ApJ, 399, 533 

Pickett, H. M., Poynter, I. R. L., Cohen, E. A., et al. 1998, JQSRT, 60, 883 

Sakamoto, K., Scoville, N. Z., Yun, M. S., et al. 1999, ApJ, 514, 68 

Sanders, D. B. & Mirabel, I. F. 1996, ARA&A, 34, 749 

Schoier, E, van der Tak, E, van Dishoeck, E., & Black, J. 2005, A&A, 432, 369 
Soifer, B. T., Neugebauer, G, Matthews, K., et al. 1999, ApJ, 513, 207 
Spaans, M. & Meijerink, R. 2007, ApJ, 664, L23 
Tilanus, R. P. J., Tacconi, L. J., Sutton, E. C, et al. 1991, ApJ, 376, 500 
Van der Tak, E, Walmsley, C, Herpin, E, et al. 2006a, A&A, 447, 101 1 
Van der Tak, F. F. S., Belloche, A., Schilke, P., et al. 2006b, A&A, 454, L99 
Van der Tak, F. F. S., Black, J. H., Schoier, F. L., et al. 2007, A&A, 468, 627 
Wilson, C. D., Booth, R. S., Olofsson, A. O. H, et al. 2007, A&A, 469, 121 



F.F.S. van der Tak et al.: Detection of extragalactic H 3 + 




Fig. 2. Calculated chemical structure of XDR (top) and PDR (bottom) models 



